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Abstract 

We undertake a comprehensive investigation of the properties of the 
sphaleron in electroweak theories with two Higgs doublets. We do this in 
as model-independent a way as possible: by exploring the physical parame- 
ter space described by the masses and mixing angles of the Higgs particles. 
If there is a large split in the masses of the neutral Higgs particles, there 
can be several sphaleron solutions, distinguished by their properties un- 
der parity and the behaviour of the Higgs field at the origin. In general, 
these solutions appear in parity conjugate pairs and are not spherically 
symmetric, although the departure from spherical symmetry is small. In- 
cluding CP violation in the Higgs potential can change the energy of the 
sphaleron by up to 14 percent for a given set of Higgs masses, with signif- 
icant implications for the baryogenesis bound on the mass of the lightest 
Higgs. 
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1 Introduction 



One of the major unsolved problems in particle cosmology is to account for the 
baryon asymmetry of the Universe. This asymmetry is usually expressed in terms 
of the parameter T], defined as the ratio between the baryon number density 
Ub and the entropy density s: 77 = %/s ~ 10~ 10 . Sakharov |lj laid down 
the framework for any explanation: the theory of baryogenesis must contain 
baryon number (B) violation; charge conjugation (C) violation; combined charge 
conjugation and parity (CP) violation; and a departure from thermal equilibrium. 
The Standard Model is naturally C and P violating, and violates CP through 
the couplings of fermionic charged currents to the (the CKM matrix). It was 
also known to violate the combination B + L (where L is lepton number) non- 
perturbatively 0, and the realisation that this rate is large at high temperature, 
and that the Standard Model could depart from equilibrium at a first order 
phase transition || led to considerable optimism that the origin of the baryon 
asymmetry could be found in known physics. 

However, the Standard Model does not have a first order phase transition for 
Higgs masses above about 75 GeV ||, and in any case is not thought to have 
enough CP violation. Current attention is focused on the Minimal Supersym- 
metric Standard Model (MSSM), where there are many sources of CP violation 
over and above the CKM matrix || [7], and the phase transition can be first 
order for Higgs masses up to 120 GeV, provided the right-handed stop is very 
light and the left-handed stop very massive || [10], pT| . 

The currently accepted picture for the way these elements fit together was 
developed by Cohen, Kaplan, and Nelson [12j (see also |L3|, [L4|, [15| for reviews). 



A first order transition proceeds by nucleation of bubbles of the new, stable, 
phase. The bubbles grow and merge until the new phase has taken over. The 
effect of CP violation in the theory is to make the fermion reflection coefficients 
off the wall chirally asymmetric, which results in a chiral asymmetry building up 
in front of the advancing wall in the fermion species which couple most strongly 
to the wall and have the largest CP violating couplings. This chiral asymmetry 
is turned into a baryon asymmetry by the action of symmetric-phase sphalerons. 

As the wall sweeps by, the rate of baryon number violation by sphalerons drops 
as the sphaleron mass increases sharply. The formation of a sphaleron is a thermal 
activation process and the rate can be estimated to go as Y s ~ exp(— E S (T)/T), 
where E S (T) is the energy of the sphaleron at temperature T. This rate must 
not be so large that the baryon asymmetry is removed behind the bubble wall by 
sphaleron processes in thermal equilibrium, and this condition can be translated 
into a lower bound on the sphaleron mass [[R], [17], [ll| 

E S (T C )/T C Z45. (1) 

Thus it is clear that any theory of baryogenesis requires a careful calculation of 
the sphaleron mass. For example, it turns out that condition (TJ) is not satisfied 
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for any value of Higgs mass in the Standard Model || . 

It has been known for a long time that spherically symmetric solutions exist 



in SU(2) gauge theory with a single fundamental Higgs [19, 2C, 21|, which is 



the bosonic sector of the Standard Model at zero Weinberg angle. However, it 
was Klinkhamer and Manton [fSJ who realised that they were unstable, with a 
single unstable mode, and that the formation and decay of a sphaleron results in 
a simultaneous change of both B and L number by Nf (the number of fermion 
families). They calculated numerically both the mass and the Chern-Simons 
number, finding the mass to be 3.7 (4.2) Mw/ctw at a Higgs mass of 72 (227) 
GeV, where aw — 9w/^ n an d Mw is the mass of the W ± particle; and the 
Chern-Simons number to be exactly 1/2. 

At Mh £ 12M\y new solutions appear |23|, 23|, which have different boundary 
conditions at the origin: the Higgs field does not vanish. These spontaneously 
violate parity and occur in P conjugate pairs with slightly lower energy than the 
original sphaleron, which correspondingly develops a second negative eigenvalue. 
These are termed deformed sphalerons or bisphalerons. 

Several authors have considered models with two Higgs doublets. Kastening, 
Peccei, and Zhang (KPZ) |^5| studied models with CP violation, but did not 
use the most general spherically symmetric ansatz, limiting themselves to a par- 
ity conserving form. Bachas, Tinyakov, and Tomaras (BTT) |26| on the other 
hand, considered a two-doublet theory with no explicit CP violation, used a C 
conserving ansatz, chose the masses of the pseudoscalar (Ma) and the charged 
Higgs (Mff±) to be zero, and chose the mixing between the two scalar Higgses to 
be zero. They found new P violating solutions, specific to multi-doublet models, 
at M H ^ 5M W , where M H is the mass of the second CP even Higgs. They did 
not calculate the Chern-Simons number, but we show that these solutions appear 
in P conjugate pairs and are in fact sphalerons, in that they have Chern-Simons 
number near 1/2, and one unstable mode. In view of the difference in behaviour 
of the two Higgs fields as the origin is approached, we call them relative winding 
(RW) sphalerons. More recently, Kleihaus |27| looked at the bisphalerons in a 
restricted two-doublet Higgs model. 

Sphalerons in the MSSM have were studied by Moreno, Oaknin, and Quiros 
(MOQ) [p8fl , who included one-loop corrections, both quantum and thermal. 
However, they again did not allow for P violating bisphalerons or RW sphalerons, 
and did not consider the effect of CP violation either, which can appear in the 
guise of complex values of the soft SUSY breaking terms in the potential. 

All of the above work was carried out at zero Weinberg angle with a spherically 
symmetric ansatz: there have been several studies of sphalerons in the Standard 
Model in the full SU(2)xU(l) theory |3(], [JT], where one is forced to adopt 
the more complicated axially symmetric ansatz: Ref. |5!| used the axially sym- 
metric ansatz in a numerical computaton, |5D[ expanded in powers of g'/g using 



a partial wave decomposition, and [|3T| estimated the energy by constructing a 
non-contractible loop in field configuration space which was sensitive to 8w- The 
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upshot of this work is that working at the physical value of the Weinberg angle 
changes the energy of the sphaleron by about 10%. It is interesting to note that 



the SU(2)xU(l) theory also contains charged sphaleron solutions [32 



Here we report on work on sphalerons in the two-doublet Higgs model (2DHM) 
in which we study the properties of sphalerons in as general a set of realistic 
models as possible, although we do use the zero Weinberg angle approximation 
and a spherically symmetric ansatz. We try to express parameter space in terms 
of physical quantities: Higgs masses and mixing angles, which helps us avoid 
regions of parameter space which have already been ruled out by LEP, or where 
the vacuum is unstable. It also means one can take into account ultraviolet 
radiative corrections by using the 1-loop corrected values for the masses and 
mixing angles. 

We are interested in the energy, the Chern-Simons number, the symmetry 
properties, and the eigenvalues of the normal modes of the various sphaleron 
solutions in the theory, as functions of the physical parameters. From the point 
of view of the computation of the rate of baryon number violation, the mass is 
certainly the most important quantity, followed by the number and magnitude of 
negative eigenvalues of the fluctuation operator in the sphaleron background: the 
largest contribution to the baryon number violation rate comes from the sphaleron 
with lowest energy and hence only one negative eigenvalue. The Chern-Simons 
number and the symmetry properties under C, P, and spatial rotations, are also 
interesting as they help classify the solutions. 

We firstly check our results against the existing literature, principally Yaffe 
^4] and BTT |26 |, and then reexamine the sphaleron in a more realistic part 



of parameter space, where Ma and Mjj± are above their experimental bounds. 
We find that in large regions of parameter space, particularly when one of the 
neutral Higgses is heavy (above about 6 Mw), the RW sphaleron is the lowest 
energy sphaleron. When there is CP violation in the Higgs sector, the would- 
be pseudoscalar Higgs can play the role of the heavy Higgs, and the other two 
Higgses can remain relatively light. The fractional energy difference between the 
RW and the ordinary (Klinkhamer-Manton) sphaleron is small, about 1% in the 
parameter ranges we explored. 

We encounter a problem with P violating sphalerons when either Ma — Mh± , 
or the amount of CP violation is non-zero: there is a departure from spherical 
symmetry in the energy density, signalling an inconsistency in the ansatz for the 
field profiles. However, the energy density in the non-spherically symmetric terms 
is small, at most about 0.2% of the dominant spherically symmetric terms, so it 
is a good approximation to ignore them. 

We also looked at the sphaleron in the restricted parameter space afforded by 
the (tree level) MSSM, confirming the results of [28| that the sphaleron energy 
depends mainly on the mass of the lightest Higgs and on tan/5, and finding no 
RW or bisphaleron solutions. 

Finally, we amplify the point made in that introducing CP violation 
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makes a significant difference to the sphaleron mass, and may significantly change 
bounds on the Higgs mass from electroweak baryogenesis. 

We do not explicitly compute quantum or thermal corrections fT8|, B3, |35 



38], |3£| as they are model-dependent. However, if particle masses are 
expressed in units of M\y, a reasonable approximation to the 1-loop sphaleron 
mass (in units of Myy/oiw) can be obtained by interpreting the masses and mixing 
angles as loop-corrected quantities evaluated at an energy scale M w [p9| . This 
approximation justifiably ignores small corrections due to radiatively induced 
operators of dimension higher than 4, but does not take into account the cubic 
term in the effective potential. This means our calculations are less accurate 
near the phase transition. However, as the error is in the Higgs potential, which 
generally contributes less than 10% to the energy, the resulting uncertainty is not 
large. 

The plan of the paper is as follows. In Section |2] we describe the bosonic 
sector of the two Higgs doublet SU(2) electroweak theory. We discuss the var- 
ious parametrizations of the scalar potential, and provide translation tables in 
Appendix We show how we use physical masses and mixing angles as inde- 
pendent parameters of the theory. Although in this approach the stability of the 
vacuum is automatic, as one chooses the masses of the physical particles to be 
real, there are still the problems of boundedness and global minimisation to be 
overcome. We solve the boundedness problem straight forwardly, but with two 
Higgs doublets, finding the global minimum of the potential is non-trivial, and 
we are forced to use numerical methods. 

In Section |] we discuss the sphaleron solutions and their symmetry proper- 
ties. In Section |] we descibe the numerical method we use to find the solutions: 
although the Newton method has been used before [26| there are some dif- 
ficulties associated with the boundary conditions that were not highlighted by 
previous authors. In Section we present our results. Section |] contains discus- 
sions and conclusions. 

Throughout this paper we use h = c = k B = 1, & metric with signature 
(+, -, -, -), and M w = 80.4 GeV. 



2 Two Higgs doublet electroweak theory 

We shall be working with an SU(2) theory with two Higgs doublets a , with 
subscript a = 1,2. Although we should strictly work with the full SU(2)xU(l) 
theory, neglecting the U(l) coupling is a reasonable approximation to make when 
studying the sphaleron. 

The relevant Lagrangian is 

C = --F" F a ^ + (DM\D^ a ) - Vfa, 2 ). (2) 
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Here, the covariant derivative D^(j) a 



dn4> a + gWj^t a (J) a with antihermitian gen- 



P : 




-> a (t, 




C : 


(p a (t,x J ) 






CP : 


4) a {t,x 3 ) 







erators t a = a a /2i. 

This Lagrangian may have discrete symmetries, including parity, charge con- 
jugation invariance, and CP [40|. These transformations are realised on the Higgs 
fields by 

(3) 
(4) 
(5) 

where 9 a are phase factors that can only be determined by reference to the com- 
plete theory. The transformations on the gauge fields are 

(6) 
(7) 
(8) 

With these transformations the only place a departure from C, P, or CP invari- 
ance can occur in Lagrangian (g) is in the Higgs potential term V(<pi, <p2)- 



P : 




-> W*(t,-x J ), 


C : 


W,(t,x j )- 


- {-ia 2 )W;{t,x?){-ia 2 )\ 


CP : 


W,(t,x j )- 


-> (-ia 2 )W**(t,-x j ){-ia 2 y 



2.1 The Higgs potential 

The most general two Higgs doublets potential has 14 real parameters, assuming 
that the energy density at the minimum is zero. We shall consider one with 
a discreet symmetry imposed on dimension four terms, <p\ — > <f>i, (j) 2 — > — 02> 
which suppresses flavour changing neutral currents [ 41 1 , and results in a potential 
with 10 real parameters. One of these parameters may be removed by a phase 
redefinition of the fields we detail in Appendix |A], and the potential may be 
written 



(Ai + A s 

+2A 3 
+A 4 
+ (A+ 
+2X2 



tt J +(A 2 + A 2 



^1^2 - Re 2 



402 



Co 



V t ) I Re(0i 



^2j " 
^1^2 

2 y 



^1^2 

2 J 

) 



) 2 + (A+-Xi)Im 2 



(9) 



This form of the potential is convenient as the vacuum configuration, which we 
take as the zero of the potential is entirely real: 



l^vac 



Vg 

s/2 




1 



(10) 
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This form also makes clear what are the sources of CP violation in the theory. 
Ignoring couplings to other fields, it can be seen that when \2 — there is a 
discrete symmetry 

<p a -> -icr 2 (f)* a , (11) 

which sends lm(0[0 2 ) — > — lm(0|0 2 )- This can be identified as charge conjugation 
invariance. Thus \2 is a C breaking parameter. In the presence of fermions, C 
and P are not separately conserved, and we generally refer to the field properties 
according to their behavior under CP, and to X2 as a CP violating parameter, 
giving rise to a mixing between the CP odd and CP even neutral Higgses. When 
one includes the other fields of the full theory one can find further sources of CP 
violation, such as the phases in the CKM matrices of the quarks and, if neutrinos 
are massive, leptons. 

In Appendix |A| we write down how the nine parameters of Eq. || relate to the 
parameters of the two more usual forms of this potential. 

It is useful to determine as many as possible of the nine parameters in the 
potential from physical ones. The physical parameters at hand are the four 
masses of the Higgs particles, the three mixing angles of the neutral Higgses, and 
the vacuum expectation value (v) of the Higgs (which is determined from M\y, 
and the SU(2) gauge coupling g). This leaves one undetermined parameter which 
may be chosen in various ways. 

In the absence of CP violation, we automatically have \2 = 0, and our input 
parameters are; v, Mh and Mh (the masses of the CP even scalars), Ma (the 
mass of the CP odd scalar), M#± (the mass of the charged scalar), <p (the mixing 
angle between the CP even scalars), tan (3 and A3, (the only parameter we choose 
by hand). This gives non-zero values for the other eight of our nine parameters. 

In the presence of CP violation our input parameters again include v, Mh, 
Mh, Ma, Mh±, <fr, and A3. However, now we also have 9c p (the mixing angle 
between the CP even and the CP odd neutral Higgs sector which is entirely 
responsible for the \2 term), and the third mixing angle ip. For a non-zero 6cp, 
tan/? is determined by the masses and mixings, and although we still denote 
the three neutral Higgs masses as Mh, Mh, and Ma we stress that they are not 
respectively CP even, CP even, and CP odd, but have some combination of 
these properties depending on the values of 9c p and 4>. 

The conversion between the parameters of Eq. ^ and these masses and mixings 
is carried out in the charged sector through 

K = —^, (12) 

v 

and in the neutral sector by writing 

v 2 X = D- 1 {^,9 C pA) M P (M h ,M H ,M A ) D(^,9 C pA), (13) 
where Mp is a diagonal mass matrix given by 

M P = Diag [M H , M^, M%\ , (14) 
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and D is the orthogonal matrix which diagonalises X. Defining rotation matrices 
in the usual way, 



R*(a) 



cos a sin a 
— sin a cos a 
1 



R y {a) 



cos a — sin a 

1 
sin a cos a 



(15) 



we can arrange for the mixing angles ip, 6cp, 4> to be the usual Euler angles, 
through 

D^, 6 CP , </>) = R z {^)Ry{9)R z {4>). (16) 

The X(i/j, 9cp, 4>, M h , M H , Ma) of Eq. [13] can be obtained as a function of the 
parameters of Eq. 0, by expanding about the vacuum state Eq. jTDl to give 



X(l,l) 


= - [4(A 1 + A 3 )cos 2 /3+(A + + xi)sin 2 /3] , 


(17) 


X(l,2) 


= X(2,l) = i(4A 3 + A+ + xi)cos/3sin/3, 


(18) 


X(l,3) 


= X(3, 1) = -X2sin/3, 


(19) 


X(2,2) 


= i [4(A 2 + A 3 )sin 2 /5 + (A + + xi)cos 2 /5] , 


(20) 


X(2,3) 


= X(3,2) = i X 2Cos/3, 


(21) 


X(3,3) 


= ^(A+-Xi)- 


(22) 



Inverting Eqs. |i~7] - |2~2l gives^ 

2v/X(l,3) 2 + X(2,3) 2 , 
arctan [X(l, 3)/X(2, 3)] , 

[X(l, 1) cos/5 - X(l, 2) sin/3 - 2A 3 cos 2(3 cos p] 



X2 



Ai 
A 2 
A + 
Xi 



1 



[X(2, 2) sin /? - X(l, 2) cos/3 + 2A 3 cos 2(5 sin /?] 



2 cos 3 /?' 
1 

2 sin 3 /?' 



-2A 3 + X(1,2)- 



+ X(3,3), 



sin (3 cos /3 
-2A 3 + X(l, 2)— i - X(3, 3), 



(23) 
(24) 

(25) 
(26) 
(27) 
(28) 



sin (3 cos /3 

where the X above are the X(if), 9cp, <t>, Mh, Mh, Ma) as given by Eq. [13|. And 
we have chosen — it < 2(3 < tc from which, depending on the sign of X(l, 2) and 
X{1, 3), we can set the sign of \2- Although it is unconventional to allow (3 to 
take negative values, it is a natural consequence of allowing the mixing angles to 
vary over their full range. 

1 We have corrected two typographical errors from j33|: a swapped cos and sin in Eq. p5| 
and Eq. E6L and a sign error in Eq. pq. 
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2.2 Boundedness and stability of the Higgs potential 



Before proceeding, we re-examine the conditions on our potential which derive 
from its boundedness and the stability of the vacuum state. For boundedness we 
need consider only the quartic terms of Eq. ^| to find the large field behaviour of 
the potential. We write our doublets as 



cos p a e 
sin p a e u 



(29) 



this will allow us to express the potential in terms of independent quantities. The 
quartic terms of Eq. can then be written as 



V = aQj + bQ 2 2 + c(rn,ri 2 )QiQ 2 , 



where 



V2 



cos pi cos p2 cos(k 2 — Ki) + sin pi sin p 2 cos(cj 2 — ui), 
cos pi cos p2 sin(/t 2 — «i) + sin pi sin p 2 sin(u; 2 — u>i), 



(30) 



(31) 
(32) 



and 



a 
b 



(33) 
(34) 
(35) 



Ai + A 3 , 
A 2 + A3, 

c(??i, 772) = 2A 3 + A 4 + (A+ - A 4 + Xi)Vi + (A+ - A 4 - xi)vi + 2 X2ViV2- 

The variables Qi, Q2, 771, and r? 2 are then independent. Furthermore, Qi and Q 2 
are by definition non-negative, and 7]i and r/ 2 are constrained to lie in the unit 
disc 

0<77i 2 +% 2 <1- (36) 

The potential can now be viewed as a quadratic form in Q\, Q 2 , in which case 
the form must be positive for all values of 771, 772 in the unit disc. If 0^^(771, 772) 
is the minimum value of 0(771,772) for all rji and 772, the condition for the form to 
be positive and the potential bounded are 



ab 



a + b > 

-.2 

in ~> 



0. 

0. 



(37) 
(38) 



On substituting the values of a, b, and c m i n into Eqs. |37] and [$8| we obtain 

Ai + A 2 + 2A 3 > 0, (39) 
4A 4 A 2 + 4(Ai + A 2 )A 3 - (4A 3 + Ac)A c > 0, (40) 

where 



A 



c 




if A+ 
otherwise 



Vxi 



xl 



> \ A 



(41) 



Eqs. |39] and ^ are the necessary and sufficient conditions for a bounded quartic 
potential. In 



33 1 we considered only Eqs. 



and [40] for the second case of Eq. 
The condition for the vacuum of Eq. [TO to be a minimum is simply 



m h > 0, m H > 0, 



> 0, m H ± > 0. 



(42) 



On substituting masses and mixings from Eqs. [12] and [13], and Eqs. p3]-p8| into 
the inequalities Eqs. |3^ and [|(] we could derive six conditions directly on masses 
and mixing angles. Vice versa, by substituting the expressions for the masses in 
to the parameters of the potential, six conditions could be obtained directly on 
the parameters of Eq. |9|. In practice, we picked masses and mixings, calculated 
the parameters of Eq. pi and then verified that Eqs. [39] and HO] held. 



2.3 Global minimisation 

While the constraints of Eq. [42] guarantee that Eq. [10] is a minimum of the 
potential, they do not guarantee that it is a global minimum. We are dealing with 
a large number of parameters, and before we proceed we need to be aware that 
for some regions of this parameter space the minimum of Eq. [10] is not a global 
minimum. We were unable to find all but the simplest analytic conditions on the 
parameters of our potential that constrained Eq. to be a global minimum. 

Our approach was perforce numerical: we ran the Maple extremisation routine 
extrema which took as input parameters the masses and mixings mentioned 
above. However, we found this extremisation routine was not fully reliable and did 
not find all the extrema. We instead adapted the code written to find sphaleron 
solutions to find extrema with constant fields, and looked for configurations with 
negative energy. In Appendix [B] we give more details of our numerical method of 
finding global minima. 



3 Sphaleron ansatz and spherical symmetry 



A sphaleron is a static, unstable solution to the field equations representing the 
highest energy field configuration in a path connecting one vacuum to another. It 
is easiest to look for spherically symmetric solutions, and so we use the spherically 
symmetric ansatz of |42|], extended to allow P, C pjfl , and CP violation |33| : 




1 



1 1 a ~a° a 

—=-A x a — 
V29 2z 



(43) 
(44) 
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Wi 



1 1 



(y/2 + P) A a, . M A 

^aij&j ~t~ v^ai %a^i) ~t~ -<*-iX a Xi 



a 
2i 



(45) 



i,,, and -F^, G Q , 



a, /3, Aq, and Ai are 



where F Q = a a + z6 a and G a = c a + ic? c 
functions of the radial co-ordinate r. 

We work in the radial gauge where A± is zero, and as we are looking for static 
solutions we set Aq to zero. We have scaled separately the Higgs and gauge parts 
of this ansatz so that the kinetic contribution to the energy is of the form \ f2i 
where f& generically denotes the fields 



Under the P, C, and CP transformations of Eqs. |3|-|8], where we have set 
= 0, the fields Ja, Aq, and A\ transform as shown in Table [Q. 



P 



C 



CP 



b a 
c a 
d a 
a 

P 
A 
A 1 



+a a 
+b a 

Cq, 

da 

—a 

+P 

-A 

-A, 



b a 
c a 
d a 
a 

P 
A 
A, 



+a a 
-b a 
+c a 

da 

+a 

+P 

+A 

+A X 



O'a 
b a 
C-a 
da 

a 

P 
A 
A x 



+a a 

-b a 

+d a 

—a 

+p 

-A 
-A, 



Table 1: P, C, and CP transformations for the fields of ansatz (431445 



On substuting ansatz Eqs. |43|445| into the Lagrangian |2] we find the static 
energy functional 

E[f A ] = ^jf J drd9d<P r 2 sin6 [K + V H ] (46) 

where r is in units of , and 

K = K + K lX3 , (47) 
V H = V + V 1 x 3 + V 2 x 3 x 3 . (48) 

K , K\, Vo, V\, and V 2 are given in Appendix [C], and £3 = 2<\rf a a <$" 2 x a j v\v 2 
is the third component of a unit radial vector. Hence this ansatz is potentially 
inconsistent if K\, V\, and V 2 are non-zero. 

If the field configuration conserves C: F a = a a and G a = c a , and we have the 
usual ansatz of [J2|. This gives Ki—0 and Vi—0, although V 2 may be non-zero 



if Ma 7^ M H ±, and the field configuration has c a 7^ 0. If the field configuration 
conserves P: G a = 0, and again all three of the dangerous terms Ki, V\, and V 2 
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vanish. In the presence of two Higgs doublets Bachas, Tinyakov, and Tomaras 
[p6[| (for RWS) and Kleihaus |27j (for bisphalerons) used a C conserving ansatz 
and worked with parameters for which Ma = Mn± = and thereby conserved 
spherical symmetry. On introducing C violating terms Kastening, Peccei, and 
Zhang |25[ used a P conserving ansatz to find the ordinary sphaleron, while 
in extending to the MSSM Moreno, Oaknin, and Quiros |28| used a C and P 
conserving ansatz for the sphaleron, and so again neither |25j nor ]28| would have 
noticed any departure from spherical symmetry. 

The functions Kq, Vo, and Vi for the C conserving ansatz, and the conditions 
on parameters and solutions which conserve exact spherical symmetry are given 
in Appendix |C| If we allow an ansatz which does not conserve P, C, or CP 
F a = a a + ib a and G a = c a + id a , and Ki, Vi, and Vi can all be non zero. K , 
Ki, Vo, Vi, and Vi for this case are also given in Appendix |C[ 

Our strategy is to assume Ja = /U( r ) and integrate over x 3 = cos 6* of Eqs. 



46-48 to give 



Mi 



w 



aw 



dr 



K + V + -V 2 



(49) 



If solutions, corresponding to extrema of Eq. E^, have field profiles for which 
Ki — 0, V\ — 0, and V 2 = 0, then the solutions are exactly spherically sym- 
metric, and the ansatz has succeeded. Otherwise, the solutions are not exactly 
spherically symmetric, with Ki, Vi, and V2 measuring the departure from spher- 
ical symmetry. We can then regard Eq. ^ as the first term in an expansion in 
spherical harmonics, and our procedure finds a good approximation to the / = 
modes provided that K\, V\, and V2 are all small in comparison to K Q and Vq. 

In our previous paper |33| we assumed spherical symmetry at the level of the 
static energy functional by imposing 



A(r)G c 



(50) 



which is too restrictive when it comes to finding C and P violating solutions in 
C violating theories. 



3.1 Properties of solutions 

We can classify solutions according to which of the symmetries C, P, and CP they 
preserve. The ordinary (Klinkhamer-Manton p2|) SU(2) sphaleron preserves 
both C, and P, and its extension to a C conserving two Higgs doublet theory 
therefore has a = 0, b a = 0, c a = 0, and d a = 0. Kunz and Brihaye |23[ and 



Yaffe [24] showed that, with one Higgs doublet, there exist P violating solutions 
at large Higgs mass with lower energy than the ordinary sphaleron, this solution 
is named the bisphaleron as it occurs in P conjugate pairs. The appearance of 
a bisphaleron solution is signalled by the ordinary sphaleron developing an extra 
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negative eigenvalue as the Higgs mass increases. In a C conserving theory these 
solutions are C invariant and have b a = and d a = 0, and are distinguished from 
the ordinary sphaleron by non-zero c a and a. To date they have been investigated 
with only Mh, Mh, and tan/3 non zero, which corresponds to M#± = Ma = in 
a C conserving theory, where they maintain spherical symmetry. However with 
M H ± ^ Ma or a non-zero 6 C p; V 2 , or K\, V\, and V 2 respectively can all be 
non-zero. Hence, departure from spherical symmetry is generic, even in the pure 
SU(2) two doublet model. 



Bachas, Tinyakov, and Tomaras |26| investigated two Higgs doublets mod- 



els and found more P violating solutions at lower Higgs masses than the bis- 
phaleron. Although again occurring in P conjugate pairs, they are distinguished 
from the bisphaleron in that their boundary conditions require more than one 
Higgs doublet: the two Higgs fields have a relative winding around the 3-sphere 
of gauge-inequivalent field values of constant |0i| and |0 2 |- Thus we refer to them 
as relative winding or RW sphalerons or RWS. If we refer just to a sphaleron, 
we shall henceforth generally mean the ordinary P and C conserving sphaleron. 
Note that RW sphalerons are spherically symmetric in C conserving theories only 
when Ma = M H ±. 

The defining characteristic of a sphaleron is that it represents the highest 
point of a minimum energy path starting and ending in the vacuum, along which 
the Chern-Simons number changes by ±1. The Chern-Simons number is defined 

as 

,2 



k 



Wtd 3 W% + -ge abc W t a W°W, 



(51) 



327T 2 



d 6 xK°, (52) 



where d^K^ = F^ u F a ' 1 ' / . Under a gauge transformation, ncs changes by an in- 
teger: hence, field configurations with integer ncs are gauge equivalent to the 
vacuum Wf = 0. One should also note that n cs is odd under CP. 

Ordinary sphalerons have half-integer Chern-Simons number ncs, which by 
choice of a suitable gauge can be taken to be precisely 1/2. However, Yaffe found 
that the bisphalerons pairs had ncs — 1/2 ±z/, where v was typically fairly small, 
and depended on the parameters in the Higgs potential. Bachas, Tinyakov, and 
Tomaras did not calculate the Chern-Simons number of their relative winding 
sphalerons pairs, but we also find them to come in pairs with ncs — 1/2 ± z/. 
That solutions which spontaneously violate CP in this way should come in such 
pairs is clear, as field configurations with ncs — 1/2 — u can be obtained from 
one with ncs = l/2 + z/bya combination of a CP and a gauge transformation. 



4 Finding solutions 
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4.1 Method 



We will be finding solutions to a static energy functional of the form 

M w 



aw 



dr £{f A 



where 



(53) 



(54) 



Here, -P(/U) is a polynomial in the 10 fields /a, which we divide into gauge fields 
f G = a,P and Higgs fields f H = a a , b a , c a , d a . 

We use a Newton method, following pi| , which is an efficient way of finding 
extrema (and not just minima). The method can be briefly characterised as 
updating the fields Ja by an amount 6 /a, given by the solution of 

ojB = - — , (55) 



Sf B Sf A ' SfA 

which we can abbreviate as S"5f = —£ . Provided £ has no zero eigenvalues, the 
equation has a unique solution, subject to boundary conditions which we detail 
below. We sometimes added a fraction of 5f which, although slower, occasionaly 
produced a more stable convergence. The procedure is started from an initial 
guess for /a, and then repeated with each improved configuration, until £ is 
small enough so that 8f~0. 

A particular advantage to using this method is that because we are calculating 
£ , it is straight forward to get the negative curvature eigenvalues, u 2 , from the 
diagonalistion of £ at each solution. To achieve this we use 



from 



S 2 E[f A 



Mi 



" Sf G " 


= u 2 


' Sf G ' 


_rdf H 




_r f H _ 



w 



aw 



dr 



SfG 

rSf H 



2 



8 fa 
r8f H 



(56) 



(57) 



where it is understood that the £ of Eqs. 

respect to f G and rfn, and not as in the Newton method of Eq 
to f G and f H . 



and |57| has been differentiated with 

I with respect 



4.2 Boundary conditions 

Next we turn our attention to boundary conditions. Before we look at specific 
conditions for different solutions, we consider the terms of Eq. |114j of Appendix 
0, (up to numerical factors) 

K° oc ^(a 2 + (3 2 -2r + (al + bl + cl + dl)(a 2 + ^ + 2) 

+V2(3(a 2 a + b 2 a -c 2 a - d 2 a ) - 2V2a(a a c a + b a d a ). (58) 
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We introduce new fields x, K-a, £ a , ^, and Q defined by 

-fi + ia = v / 2xexp(z^), (59) 

V 

&a + ica — 2— /C a exp(z9 Q ), (60) 

v 

V 

b a + id a = 2—£ a exp(iQ a ), (61) 

(62) 

and rewrite (1581) as 



K oc ^(X 2 -1) 2 

+ (2 X 2 + 2) [cos 2 /3(/C 2 + £ 2 ) + sin 2 /3(/C 2 + £ 2 ) 
-4x cos 2 /?(/C 2 + £ 2 )Re[exp(-i^ + 220!)] 
-4xsin 2 /?(/C 2 + £ 2 )Re[exp(-^ + z20 2 )]. (63) 

We have a boundary condition from the finiteness of the energy density, due to 
the first term in Eq. ^ which can be expressed as 

X 2 -> 1 as r -> 0. (64) 

From the finiteness of the gauge current density (which is proportional to the 
second, third, and fourth terms in Eq. ^) and using Eq. 0, we also have 



(ICj + £ 2 )Re[exp(-^ + i26i)] -> /C 2 + Cj 
(/C 2 + £ 2 )Re[exp(-^ + z26 2 )] -> /C 2 + £ 2 



as r — > 0. (65) 



To satisfy Eq. B5| we require 



/C 2 + £ 2 ^0 \ 0x ^^/2 + rixTr 



dther K 2 + £ 2 - j ° r 2 - ¥/2 + n 2 7r > aS r -> °' ( 66 ) 



where ni, n 2 € Z. Eq. |66| can be rewritten as 

either /C 2 + £ 2 ^0 

or 0i-02 — > (ni - n 2 )7r 



as r — > 0. (67) 



Eqs. and are then our boundary conditions as r — > 0. The boundary 
conditions as r — >• oo can be obtained from finiteness of K , (Eq. |107|) and of V 
(Eq. 03)- 

The ordinary sphaleron satisfies Eq. |67| by having 

(/C 2 +£ 2 )| r=0 = 0. (68) 
The full set of boundary conditions for the sphaleron are given in Table |2|. 
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Bisphaleron pairs have different boundary conditions. To satisfy Eq. |J7|, where 
5 is a small positive angle, they have 



20! 



r=0 



26 



2\r=0 



r=0 



26 



-TV ±5. 



(69) 



The boundary conditions on the Ja of these solutions are given in Table 
Relative winding sphalerons pairs satisfy Eq. ^ through 



2(6! 



7T 



r=0 



26 



2\r=0 



r=0 



-TV 



±5. 



(70) 



From Eq. |67] we see that since n\ = n 2 for bisphalerons while n\ = ^2 + 1 for RWS, 
RWS unlike bisphalerons can only occur in multi-doublet theories. The integers 
ni and n 2 represent the winding numbers of the Higgs fields around the 3-spheres 
of constant |</>i| and \4>2\i with only their difference having any gauge- invariant 
meaning. The RWS boundary conditions are given in Table 01. 







a 



(3 -> y/2 







c n 



d c 







a 



/? -y/2 



ai 
a 2 



2cos/3 
2 sin/? 



b 2 








Cl 

c 2 








^2 



Table 2: Boundary conditions for the ordinary C, and P conserving sphaleron. 



a — > v^2 sin 26 ai 
/3-> -v^ cos 26 b x 

Cl 



2/Ci cos (3 cos 6 
2£x cos (3 cos 6 
2/Ci cos (3 sin 6 
> 2£i cos (3 sin 6 



a>2 
b 2 
c 2 

^2 



• 2/C2 sin /3 cos 6 
■ 2£ 2 sin (3 cos 6 

• 2/C 2 sin /? sin 6 
> 2£ 2 sin /3 sin 6 



Table 3: Boundary conditions at the origin for the (P violating) bisphaleron. 
The boundary conditions at infinity are the same as for the sphaleron, Table 13. 




-\/2cos^ 



a% — > 2/Cx cos (3 cos 61 a 2 

bi — ► 2£ x cos /3 cos 61 6 2 

ci — > 2/Ci cos /3 sin 61 c 2 

di — * 2£i cos/? sin 61 d 2 



2/C 2 sin (3 cos 6 2 
2£ 2 sin (3 cos 6 2 
2/C2 sin (3 sin 62 
2£ 2 sin (3 sin 62 



Table 4: Boundary conditions at the origin for the (P violating) RWS. The 
boundary conditions at infinity are the same as for the sphaleron, Table 0. 



4.3 Numerical performance 

The details of the implementation of the algorithm and the boundary counditions 
are relegated to Appendix [D]. We checked the accuracy of our code by evaluating 
the energy, negative curvature eigenvalues and Chern-Simons number for some 
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-u 2 


,2 
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ncs 
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&RWS 


,2 

-u x 







o.oy i 


















6 


4.531 


6.217 


0.279 












4.528 


5.171 


7 


4.609 


7.171 


1.225 












4.587 


4.147 


10 


4.778 


11.22 


5.962 












4.668 


3.090 


13 


4.888 


17.70 


13.27 


0.316 


4.886 


11.86 


6.546 


0.454 


4.700 


2.773 


15 


4.942 


23.49 


19.49 


0.926 


4.930 


8.447 


2.349 


0.428 


4.711 


2.670 


30 


5.147 


101.4 


98.55 


3.212 


5.031 


5.207 




0.387 


4.734 


2.451 


50 


5.243 


292.7 


290.1 


4.734 


5.052 


4.874 




0.380 


4.738 


2.403 



Table 5: Energy (Mw/ctw), negative eigenvalues (Myy), and Chern-Simons 
number for m = Mh/M^/ = Mh/M w and tan (3 = 1, for some of the same param- 
eters as |24]] and ||26|| . The solution with energy was reached by perturbing 



-u 



3 j 



the ordinary sphaleron in the direction of the eigenvector with eigenvalue 
and the solution with energy Erws was reached by a perturbation with eigen- 
value — u>2- If we refer to Fig. 2 of |26[ we see that the bisphaleron branch itself 



bifurcates at the point where it no longer has two negative eigenvalues, and we 
note as a point of interest that the eigenvector with eigenvalue — to\ takes us to 
the solution with lowest energy and not the S\ of [p6fl . The n C s of the RWS for 
equal CP even Higgses, and tan (3 = 1 with all other parameters zero is 1/2, this 



is not the case generally. The agreement with |24j] and |26j is excellent. 



of the same parameters as Yaffe []24J and Bachas, Tinyakov, and Tomaras [26|| . 
and found good agreement. These can be seen in Table |5|. 

The numerical scheme worked excellently, with typical convergence after five 
to fifteen iterations of 1 x 10~ 13 in the sum of absolute change in all fields at 
all points. The few problems we did encounter were: (1) sometimes the initial 
configuration for a RW sphaleron was so close to the sphaleron that the Newton 
extremisation found the original sphaleron, particularly at points in parameter 
space near the bifurcation point, and (2) the Newton extremisation sometimes 
found the vacuum from the inital configuration for a RW sphaleron . The first was 
solved by using a higher mass RW sphaleron as initial conditions for minimisation, 
and the second problem by updating each minimisation not with 5f a but with a 
fraction of it. 

We ran simultaneously two codes. One with the C conserving ansatz, and the 
other with the C and P violating ansatz. In the absence of C violation the two 
codes were identical. With 101 point instead of 51, the difference in energy, ncs, 
and eigenvalues was at most of order 0.5 % of the value with 51 points. 



5 Results 
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5.1 No CP violation, M A = M H ± = 

In order to compare with previous work, we firstly examine the unrealistic limit 
of Ma = M H ± = 0, with no explicit CP violation in the potential. We set the 
parameters A3 = and tan j3 = 6, and scanned through and Mh between 
and 800 GeV. 



Energy of sphaleron and RWS 




M h (Gev) 



Figure 1: Contours in M h , M H space of the energy of the sphaleron (dashes), and 
of the RWS (solid), in units of M w /aw Below the dotted line the sphaleron is the 
only solution. Above the dotted line, both solutions exist. The input parameters 
are tan (3 = 6 with all other parameters zero. 



Figs. [T]-[^ show contours in the Mh and Mh plane. The contours are respec- 
tively of energy (Fig. [I]), most negative eigenvalue and second most negative 
eigenvalue (Fig. ||), and n C s (Fig. |3]) of the sphaleron and relative winding 
sphaleron. When we show equal contours of both solutions the sphalerons are 
shown as dashes, and the RWS as solid. Below the black horizontal dotted line, 
shown on all four contour plots, only the sphaleron solution exists, above the black 
dotted line both solutions exist. The sphaleron never develops a third negative 
eigenvalue, nor the RWS a second negative eigenvalue. The solutions maintained 
exact spherical symmetry: V% was zero throughout; this was expected as both 
6cp = 0, and Ma = Mh± = 0. These contours are from the same potential as 



used by BTT |26] and contain some of the parameter space they scanned. Where 
we overlap we agree with their results, and we confirm their observation that the 
second negative eigenvalue appears when one of the Higgs has a largeish mass, 
(Mh ~ 5Myy). For low values of this heavier mass the lighter Higgs needs to 
be as light as possible; i.e. for the existence of relative winding sphalerons it is 
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Figure 2: Contours in M h , M H space of the eigenvalue in units of M^. The top 
figure shows the most negative eigenvalule of the sphaleron (dashes), and of the 
RWS (solid). The bottom figure shows the second most negative eigenvalue of 
the sphaleron. Below the dotted line the sphaleron is the only solution. Above 
the dotted line, both solutions exist. The input parameters are tan/3 = 6 with 
all other parameters zero. 
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preferable to have the two Higgs masses, Mh and Mh, well separated. 

Fig. p] shows both the energy of the sphaleron and the energy of the RWS, there 
is almost no difference between their energies, and the energy depends mainly on 
the mass of the lighter Higgs. Figure [| shows the most negative eigenvalue of both 
the sphaleron and RWS, and we see that there is a large difference between the 
values of negative eigenvalues for the different solutions; the negative eigenvalue of 
the sphaleron can be double that for the relative winding sphaleron for the same 
point in parameter space. Fig. |2| also shows the second negative eigenvalue of 
the sphaleron. The second most negative eigenvalue belongs to the perturbation 
which leads to the RW sphaleron in configuration space. 



Chern-Simons number of RWS 




200 ■ 
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°0 100 200 300 400 500 600 700 800 
M h (Gev) 

Figure 3: Contours in Mh, Mh space of the Chern-Simons number of the RWS. 
Below the dotted line only the sphaleron solution exists, with rics — 1/2. The 
input parameters are tan (3 = 6 with all other parameters zero. 

Looking at Fig. ^| we see that the Chern-Simons number of the RW sphaleron 
is generally not a half. There is a line in the contour space where rics =1/2. This 
occurs, for tan/? = 1, along the line of Mh = Mh, and shifts in the contour plane 
for different values of tan/?. We have only shown here solutions with rics < 1/2. 
Each of these solutions with rics < 1/2 has a P conjugate partner, with Chern- 
Simons > 1/2, such that rics + n^g = 1. 

5.2 No CP violation, M A = 3M W , M H ± = 2M W 

Figs. |]-^ show contours in M h , M H space of energy (Fig. most negative 
eigenvalue of the sphaleron and RWS, (Fig. [| top), and second most negative 
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eigenvalue of the sphaleron (Fig. |5| bottom), and rics (Fig. §]) °f the sphaleron 
and the relative winding sphaleron. Again when both solutions are shown the 
sphaleron is dashes, and the RWS solid. 
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Figure 4: Contours in Mh, Mh space of the energy of the sphaleron (dashes), and 
of the RWS (solid), in units of Mw/&w- Below the dotted line the sphaleron is 
the only solution, while above, both solutions exist. For the dotted area the 
potential is unbounded. The input parameters are tan (3 = 6, Ma = 241 GeV, 
M H ± = 161 GeV, and A 3 = -0.05. 

For these figures we took Ma = 241 GeV, M#± = 161 GeV, again with no 
explicit CP violation. We set the parameters A3 = —0.05, and tan/3 = 6, and 
scanned through Mh and Mh between and 800 GeV, with 20 GeV increments. 
Again below the black dotted line, shown on all four contour plots, only the 
sphaleron solution exits, while above both solutions exist. We see that the RW 
sphaleron solutions still persist for a large region of the parameter space. The 
dotted region at low Mh was unbounded according to Eqs. |39| and [4(]. These 
solutions did not maintain exact spherical symmetry corresponding to V% = 0, 
but the maximum value of energy due to the V2 term was 0.6% of the energy due 
to V Q . 

The solutions have the same general features as those at zero Ma and M H ±'- 
the RW sphaleron appears at widely separated Mh and Mh- While the energies 
of the two solutions in Fig. [| are almost indistinguishable, the most negative 
eigenvalue (Fig. |5| top), of the sphaleron can be double that of the RW sphaleron. 
We show the value of the second most negative eigenvalue of the sphaleron in 
Fig. H (bottom). The sphaleron never developed a third negative eigenvalue, nor 
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Most negative eigenvalue of sphaleron and RWS 
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Figure 5: Contours in M h , M H space of eigenvalues in units of M^. The top 
figure shows the most negative eigenvalue of the sphaleron (dashes), and of the 
RW sphaleron (solid). The bottom figure shows the second most negative eigen- 
value of the sphaleron. Below the dotted line the sphaleron is the only solution. 
Above the dotted line, both solutions exist. For the dotted region the potential 
is unbounded. The input parameters are tan/3 = 6, Ma = 241 GeV, M#± = 161 
GeV, and A 3 = -0.05. 
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Chern-Simons number of RWS 
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Figure 6: Contours in Mh, Mh space of the Chern-Simons number of the RW 
sphaleron. Below the dotted line only the sphaleron solution exists, with ncs — 
0.5. For the dotted region the potential is unbounded. The input parameters are 
tan/5 = 6, M A = 241 GeV, M H ± = 161 GeV, and A 3 = -0.05. 

the RW sphaleron a second negative eigenvalue. In Fig. ^| we show the Chern- 
Simons number of the RW sphaleron, and again for every solution shown with 
ncs = 1/2 — v there is a P conjugate solution with riQg = 1/2 + v. 

5.3 CP violation, M A = 8M W , M H ± = 2M W 

Figs. [7|-|9] show contours in Mh, Mh space of energy and second negative eigen- 
value (Fig. [7]), most negative eigenvalue (Figs. and Chern-Simons number (Fig. 
H) of the sphaleron and relative winding sphaleron. Sphaleron contours are shown 
as dashed lines and RW sphaleron contours as solid when present on the same 
graph. 

For these figures we took Ma = 643 GeV, Mh± = 161 GeV, this time with 
CP violation: 9qp = 0.49rr. The remaining parameters were = O.lrr, '0=0.0, 
and A 3 =3.0, giving tan/3=3.1. We scanned through M h between and 400 GeV, 
and M H between and 800 GeV, with 20 GeV increments. The dotted region at 
low M h was unbounded according to Eqs. [39] and J5| and for the white out area, 
surrounded by the solid black line, the minimum of Eq. |H] was not the global 
minimum. 

As with the previous contour plots, a large region of parameter space con- 
tained relative winding sphalerons. For these input parameters, though, due to 
the large CP violating mixing angle, the role of the large Higgs mass Mh is taken 
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Energy of sphaleron and RWS 
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Figure 7: Top: contours in M h , M H space of energy in units of M w /a\y of 
the sphaleron (dashes), and of the RW sphaleron (solid). Bottom: contours in 
M h , M H space of second negative eigenvalue (Myy) of the sphaleron. Above 
the dotted line the sphaleron is the only solution, while below both solutions 
exist. For the blank area Eq. ID is not the global minimum. For the dotted area 
the potential is unbounded. The input parameters are Oqp = 0.497T, = 0.l7r, 
ijj = 0.0, M A = 643 GeV, M H ± = 161 GeV, and A 3 =3.0. tan (3=3.1. 
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Most negative eigenvalue of sphaleron 
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Figure 8: Contours in M h , M H space of the most negative eigenvalue (M^) of 
the sphaleron (top) and of the relative winding sphaleron (bottom). Above the 
dotted line the sphaleron is the only solution, while below both solutions exist. 
For the blank area Eq. |H] is not the global minimum. For the dotted area the 
potential is unbounded. The input parameters are 9cp=0.4Qtc, 0=0. In, -0=0.0, 
M A = 643 GeV, M H ± = 161 GeV, and A 3 =3.0. tan 0=3.1. 
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Chern-Simons number of RWS 
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Figure 9: Contours in Mh, Mh space of the Chern-Simons number of the RWS. 
Above the dotted line only the sphaleron solution exists, with rics = 1/2- For the 
blank area Eq. [L(] is not the global minimum. For the dotted area the potential 
is unbounded. The input parameters are 9cp=0A9n, 0=0. Iti, ip=0.0, Ma = 643 
GeV, M H ± = 161 GeV, and A 3 =3.0. tan 0=3.1. 



on by Ma- Since, from previous contour plots, the relative winding sphaleron so- 
lution prefers regions of parameter space where there is a large separation in 
values of the heaviest (in this case the Ma) and the lightest (in this case Mh, and 
Mh) Higgs masses, the relative winding sphaleron solultions exist for the lower 
part of the contour plot, and not the upper part. Referring to Figs. above 
the black dotted line the sphaleron is the only solution, while below the black 
dotted line both the sphaleron and the relative winding sphaleron exist, this is 
opposite to the behaviour in the absence of CP violation. 

From Fig. [F] (top) the energy of the two solutions is as before almost the 
same. The second negative eigenvalue of the sphaleron is shown in the lower half 
of Fig. [7|. The sphaleron does not develop a third negative eigenvalue, nor the RW 
sphaleron a second negative eigenvalue. We show the most negative eigenvalue of 
the sphaleron and the RW sphaleron (Fig. ||) on separate graphs, and again their 
respective negative eigenvalues can be very different at the same point in the 
contour plane. We then show the Chern-Simons numbers for the RW sphaleron 
in Fig. Note that we only show solutions with rics < again, there are 
parity conjugate partners to each of these RW sphalerons, and the rics of the 
RW sphaleron and of its parity partner add up to one. 

There is no breaking in the degeneracy of the relative winding sphaleron pairs 
in energy, eigenvalues, or absolute difference from 1/2 of Chern-Simons number, 
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due to the presence of CP violation. The solutions are not exactly spherically 
symmetric, and have non zero values for all three of K\, V\, and V^- The values 
of K\, Vi, and V2 as a percentage of the Higgs potential energy are each never 
more than 0.5 %. 



5.4 MSSM parameter space 

Next we scan through tree level MSSM parameter space. Fig. ^ shows the scan 



in Ma, tan/3 space. Fig. |TT] shows the scan in Mh, Mh space. We plot contours 
of energy (top) and negative eigenvalue (bottom) for each of these scans. 

For the range of parameters we show the sphaleron did not develop a second 
negative eigenvalue. There was no departure from spherical symmetry, as only 
the a a field of the Higgs ansatz and the j3 field of the gauge ansatz were ever 
non- zero. From these four contours (Figs. [TO and O) we agree with the general 



result of |2g] that the energy of the sphaleron is sensitive to mainly Mh and 
tan/?, although their results should be more accurate as they included 1-loop 
radiative corrections. There were no relative winding sphalerons for the range of 
parameters explored. 

5.5 Sphaleron energy and CP violation 

We recall that a CP violating mixing angle can have a large effect on the prop- 
erties of the sphaleron. Here (Fig. [12]) we scan through M h , Q C p space and show 



the energy of the sphaleron and the negative eigenvalue of the sphaleron for in- 
put parmaters 0=O.125vr, ^=0.0, M H = 110 GeV, M A = 500 GeV, M H ± = 500 
GeV, and A3=0.0, these give tan/?=2.4. For the dotted region at low Mh the 
potential was unbounded, and for the blank region, bordered by the solid black 
line, the minimum of Eq. [l(] was not the global minimum of the static energy 
functional. For this region of parameter space the sphaleron never developed a 
second negative curvature eigenvalue. 

The energy of the sphaleron (Fig. [12|: top) is dependent upon the value of 
the CP violating mixing angle, and changes by about fourteen percent as the 
mixing angle varies between its minimum and its maximum. The energy is, in 
the presence of CP violation, still sensitive to the lightest Higgs mass. 

The negative eigenvalue (Fig. [12]: bottom) also has this strong dependence 
on the CP violating mixing angle, with an increase of over fifty percent as the 
mixing angle varies. Also the dependence on Mh, although not as dramatic as 
the effect of CP violation, is still present. 

5.6 Field profiles 
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Figure 10: Contours in M4, tan/3 space of the sphaleron for tree level MSSM 
parameters. The top figure shows energy (M w /aw) of the sphaleron. The bottom 
figure shows negative curvature eignevalue (Myy) of the sphaleron. 
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Curvature eigenvalue of sphaleron for tree level MSSM parameters 
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Figure 11: Contours in M h , M H space of the energy (M w /otw) of the sphaleron 
for tree level MSSM parameters. The top figure shows energy (M w /a\y) of the 
sphaleron. The bottom figure shows negative curvature eignevalue (Myy) of the 
sphaleron. 
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Figure 12: Top: contours in M h , 6 C p space of the energy (M w /a w ) of sphaleron, 
and bottom: of the negative eigenvalue of the sphaleron (M^). For this region 
of parameter space the sphaleron is the only solution. For the blank region Eq. 
[P0| is not the global minimum. For the dotted region the potential is unbounded. 
The input parameters are 0=O.125vr, ^=0.0, M H = 110 GeV, M A = 500 GeV, 
M H ± = 500 GeV, and A 3 =0.0. tan (3=2 A. 
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5.6.1 Sphaleron and RW sphaleron 

Next we show the field profiles for the sphaleron, relative winding sphaleron, and 
conjugate relative winding sphaleron for a point in the contour plot of Section 
|5.3| corresponding to a CP violating theory with M A = 8Mw, Mjj± = 2Mw, 
Mh = 1.25Mjy, and Mh = 1.5Mw We recall that the mixing angles were 
C p=0A9ir, 0=0. l7r, ^=0.0, and the coupling A 3 = 3.0. 

Before we proceed we check whether this point in parameter space is phe- 
nomenologically viable at zero temperature, as Mh = 1.25Mw is ruled out if the 
hZZ coupling is too large. We calculate the couplings gtzz, gnzz, and g A zz 
according to p3] using the values of input parameters used in Figs. P^P^, and 



compare them with the latest particle data |f44 
Using 



g hZZ = D[l,l] cos/? + P[2,1] sin/3 (71) 
9hzz = L>[l,2]cos/3 + P>[2,2]sin/3 (72) 
g AZZ = P[l,3]cos/3 + P[2,3]sin/3 (73) 



where D is given by Eq. [TBI, we obtain, for the paramaters of figures |Hi|-[TE 



glzz = 0.081 (74) 
g 2 HZZ = 0.824 (75) 
g 2 AZZ = 0.095 (76) 

which for masses M h = 101 GeV, Mh = 121 GeV, and M A = 643 GeV are with in 
experimental bounds. Although we have labelled the Higgses with subscripts h, 
H, and A; because of the values of the mixings </> = O.Itt, 9qp = 0.497T, ip = 0.0, 
while the particle with subscript h is CP even, those with subscript H, and A 
are a mix of CP even and CP odd. 

We then plot the energy density of the two types of solution, and the values of 
Ki, Vi, and Vi as a function of the rescaled radial co-ordinate for the sphaleron, 
RW sphaleron, and conjugate RW sphaleron. We recall that the departure of Ki, 
V\, and V*j from zero signals the breakdown of the spherically symmetric ansatz, 
and their size relative to the total energy density indicates the seriousness of the 
breakdown. 

It is convenient to plot the field values rescaled according to 

f "-i f f' (77) 

as then the asymptotic values are either or ±1. 



The ordinary sphaleron field profiles are plotted in Fig. as a function of 
the rescaled radial points. The solution has non zero values of a a , b a , and (3 as 
expected for a field configuration that preserves P but violates C, due to the 
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Figure 13: The sphaleron field profiles (top), and the profiles for b\ and 62 i n 
more detail (bottom). c a = d a = a = 0. This configuration has energy=4.053 
M w /aw, ncs=^l^i an d two negative curvature eigenvalues — 8.696M^, and 
-1.754M^. Input parameters are: 8 C p=0A9ir, 0=O.lvr, ^=0.0, M h = 101 GeV, 
M H = 121 GeV, M A = 643 GeV, M H ± = 161 GeV, and A 3 =3.0. These give 
tan/3=3.1, Ai = 26.29, A 2 = -2.59, A+ = 0.91, A 4 = 0.85, xi = 0.42, and 
X2 = 0.41. 
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Figure 14: The RW sphaleron field profiles (top) and the profiles for b\, b 2 , 
c 2 , di, d 2 , and a in more detail (bottom). This configuration has energy=4.047 
M w /aw, ^cs=0.478, and one negative curvature eigenvalue — 3.637M^. Input 
parameters are: cp =O.49tt, 0=0. lvr, ^=0.0, M h = 101 GeV, M H = 121 GeV, 
M A = 643 GeV, M H ± = 161 GeV, and A 3 =3.0. These give tan/?=3.1, Ai = 26.29, 
A 2 = -2.59, A+ = 0.91, A 4 = 0.85, xi = 0-42, and X2 = 0.41. 
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Figure 15: The conjugate RW sphaleron field profiles (top), and the pro- 
files for bx, b 2 , C2, d\, d 2 , and a in more detail (bottom). This configuration 
has energy=4.047 M w /a\y, ^cs=0.522, and one negative curvature eigenvalue 
-3.637M^. Input parameters are: CP =O.497r, 0=O.lvr, ^=0.0, M h = 101 GeV, 
M H = 121 GeV, M A = 643 GeV, M H ± = 161 GeV, and A 3 =3.0. These give 
tan 0=3.1, Ai = 26.29, A 2 = -2.59, A+ = 0.91, A 4 = 0.85, xi = 0.42, and 
X2 = 0.41. 
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Total and Higgs energy density for sphaleron and RWS 
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Figure 16: The top of the figure shows the total and the Higgs potential con- 
tribution to energy density in units of M^/aw for the sphaleron (solid) and the 
RWS (dashes). The bottom figure shows K\, ]/%, and V 2 for the RWS (solid) 
and its conjugate (dashes) in the same units. Both K\ and V\ are equal to their 
values for conjugate solutions, but have opposite sign. V2 is equal to its value 
for the conjugate solution. Input parameters are: 9cp=0A9ti, 0=0. lit, ip=0.0, 
M h = 101 GeV, M H = 121 GeV, M A = 643 GeV, M H ± = 161 GeV, and A 3 =3.0. 
These give tan/3=3.1, Aj = 26.29, A 2 = -2.59, A + = 0.91, A 4 = 0.85, Xi = 0.42, 
and X2 = 0.41. 



35 



presence of a C violating parameter in the potential. The sphaleron has Chern- 
Simons number 1/2, two negative eigenvalues (-8.696 M w , and -1.754 M w ), and 
has energy 4.053Mw/cew- 

The relative winding sphaleron field configurations, shown in Fig. [14], have 
non zero values for all fields. The solution violates P spontaneously and C ex- 
plicitly, and violates the combination CP. It has one negative eigenvalue (-3.637 
M w ), energy less than its ordinary sphaleron (4.047 M w / aw), and Chern-Simons 
number 0.478. Its parity congugate partner, shown in Figure [15], has field profiles 
identical to a P transformation of the RWS: that IS Cq/ ^ - d a — > — d a , and 
a — > —a, with all other fields remaining unchanged. The solution has identical 
energy, and eigenvalue to its P conjugate solution, and its Chern-Simons number 
is 0.522. 

Next we show (Fig. |16|: top) the energy density of the sphaleron, and the 
RW sphaleron, and in detail (Fig. [Ll]: bottom) the values of K x , V%, and V 2 for 
the RWS in units of energy density. K\ and V\ are equal in value, but opposite 
in sign for the conjugate pair, V2 is equal in value and equal in sign. These 
deviations from spherical symmetry are of order one part in 10 3 for these values 
of parameters. 



5.6.2 Bisphaleron 

For completeness we detail the bisphaleron fields profiles for non zero Ma and 
Mh± , and show their departure from spherical symmetry. Figs. [17] and [l^ con- 
cern this bisphaleron. We have chosen masses which are perhaps unrealistically 
large, in order to reach the part of parameter space where the bisphaleron exists: 
tan/3=6.0, M h =15.0M w , M H =17.0M W , M A =2.0M W , M H ±=3.0M W and A 3 =- 
0.1, with no CP violation. For these input parameters Ai = 567.6, A2 = 12.4, 
A+ = 0.627, A 4 = 1.923, \i = -0.227, and \2 = 0.0. 

The energy density and departure from spherical symmetry are shown in 



Figure |17[ The CP invariance means that b a = d a = 0, and hence K\ and V\ 
vanish. The departure from spherical symmetry is entirely in the V2 term shown 
in units of energy density (M w /a>w) in the lower half of Fig. [17]. The departure 
from spherical symmetry is of order 1 part in 10 



1 



The configuration in Fig. [L^ has energy=4.932 M w /aw, n cs =0.569, it has 
two negative curvature eigenvalues -11.915 M w , and -6.788 M w . Its associated 
sphaleron has energy=4.943 Mw/&w with ncs=l/2, and three negative curva- 
ture eigenvalues -23.823M^, -13.249M^, and -0.933M^. Its conjugate bis- 
phaleron has identical energy, and negative curvature eigenvalues, but ncs=0.431; 
so again the ties of the bisphaleron and its conjugate add to one. 
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Total and Higgs energy density for sphaleron and bisphaleron 
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Figure 17: The top figure shows total and Higgs potential contribution to energy 
density (M^/c^w) f° r the sphaleron (solid) and the bisphaleron (dashes). The 
bottom figure shows V*j for the bisphaleron solution and its conjugate. V% for 
both the bisphaleron and conjugate solution are equal. Input parameters are: 
tan/3=6.0, 6 C p=0-0, 0=0.0, ^=0.0, M h =15.0M w , M H =17.0M W , M A =2.0M W , 
M H ±=3.0M W , and A 3 =-0.1. 
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Bisphaleron field profiles 
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Figure 18: The bisphaleron field profiles for tan/3=6.0, Mh=l5.0Mw, 
M H =17.0M W , M A =2.0M W , M H ± = 3.0M W and A 3 =-0.1. It has energy=4.932 
Mw/cew, ncs=0.569, two negative curvature eigenvalues — 11.915Mjy, and 
— 6.788M^. Its conjugate partner is the identical solution under P conjugation 
(a — > —a), and has n C s=0A31. 

6 Conclusions 



In this paper we have made a thorough study of the properties of sphalerons in 
two Higgs doublet SU(2) gauge theories. Using a spherically symmetric approx- 
imation, we have performed scans in the physical parameter space defined by 
the masses and mixing angles of the Higgs particles, recording the energy, lowest 
eigenvalues, and the Chern-Simons number, with results recorded in Figs. |I|-|T2" 



We have also shown the profiles of the fields of our ansatz for selected solutions 

in Figs, ra-irs. 



We can draw a number of broad conclusions from these results. Firstly, for a 
wide range of parameters, the minimum energy sphaleron is not the natural gen- 
eralisation of the Klinkhamer-Manton sphaleron [^] with vanishing Higgs fields 
at the origin, but a parity violating pair of relative winding (RW) sphalerons, 



first identified by Bachas, Tinyakov, and Tomaras p6|. These are related to the 



bisphalerons or deformed sphalerons found in one doublet models by Yaffe [p4 
and Kunz and Brihaye 



but are specific to two Higgs doublet models. This 
pair was always degenerate in energy, as is to be expected from a parity conserv- 
ing Lagarangian. This degeneracy is lifted when Standard Model fermions are 



included f45fl . 

The favoured regions of parameter space for RW sphalerons to exist are those 
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where there is a large difference in the masses of the neutral Higgses. The mass 
of the heavier Higgs can be as low as bMyy. Bisphalerons appear at yet higher 
heavy Higgs masses, but were always more massive than the RW sphalerons in 
the parameter space we explored. 

The appearance of extra sphaleron solutions is signalled by the ordinary 
sphaleron developing another negative eigenvalue: thus where the RW sphaleron 
exists the ordinary sphaleron has two negative eigenvalues, and three where the 
bisphaleron exists also. The lowest energy sphaleron must have exactly one neg- 
ative eigenvalue. The numerically calculated eigenvalues of a solution not only 
aid its identification, but are important for accurate calculation of the baryon 
number violation rate: if the negative eigenvalue of the lowest energy sphaleron 
solution is ou^, then the rate is proportional to |u;_| [[34]. The difference between 
the most negative eigenvalue of the sphaleron and the negative eigenvalue of the 
RW sphaleron could be well over a factor of two. 

The most important quantity for the calculation of the B violation rate is 
normally the sphaleron energy. There is however very little difference in the 
energies of the ordinary and RW sphaleron: typically less than 1% in the range 
of parameters we surveyed. Thus the main contribution to the error in the rate 
from using the ordinary sphaleron comes from the negative eigenvalue. One 
must not only use the correct eigenvalue but also include a factor of two in the 
RW sphaleron rate, one for each of the two degenerate parity conjugate solutions. 
However, this leads only to logarithmic corrections to the sphaleron energy bound 

The most important parameter for the sphaleron energy was found to be 
the mass of the lightest Higgs, in accordance with previous studies. However, 
we were able to extend our work on the dependence of the energy on the CP 
violating mixing angle 8c p |33| to show that there was an strong dependence on 
this quantity as well, with the sphaleron energy varying by ~ 15 % as Qqp was 
adjusted through its allowed range. We note as well that we were unable to find 
a region of parameter space for which RW sphalerons existed over a wide range 
of 9 C p, for which the potential was bounded, and for which Eq. |H] was the global 
minimum. 

Although we used a spherically symmetric ansatz, we found that two Higgs 
doublet sphalerons are generically not spherically symmetric. This means that 
our results are approximate: however, the departure from spherical symmetry, 
as measured by the relative size of the symmetry violating terms in the static 
energy functional, was less than 0.2%, and so this is not a serious problem for the 
accuracy of our results. A larger correction is to be expected when one considers 
the full SU(2)xU(l) theory at non-zero 8w, for which one also has to abandon 
the spherically symmetric ansatz and resort to an axially symmetric one instead 

Another source of error is the neglect of radiative and thermal corrections. Ide- 
ally one should work out the determinants of fluctuation matrices [B51 BBI, 157], BH]. 
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One can also find solutions using the 1-loop finite temperature effective potential 



2q| . This is an implicit gradient expansion, neglecting finite temperature correc- 



tions to gradient terms, which turn out to be small [39]. Such computations are 
model-dependent: one first computes radiatively corrected couplings in the static 
energy functional, and then the sphaleron energy. Our approach decouples the 
computation of the radiative corrections, for we can take masses and angles to 
be their 1-loop corrected values. Although this neglects cubic terms and terms of 
dimension higher than 4 in the potential, it is an easy way of improving on the 
tree-level calculation, without sacrificing too much accuracy, as the contribution 
to the energy from the Higgs potential can be seen from Figs. [16] and |D] to be 
small. 

Despite these sources of error, we can conclude the calculations of the sphaleron 
energy in CP conserving models cannot safely be applied to CP violating elec- 
troweak theories, and that the sphaleron bound on the mass of the lightest Higgs 
in CP violating theories requires further investigation. 

We wish to thank Mikko Laine and Neil McNair for helpful discussions. This 
work was conducted on the SGI Origin platform using COSMOS Consortium fa- 
cilities, funded by HEFCE, PPARC and SGI. We acknowledge computing support 
from the Sussex High Performance Computing Initiative. 
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A Parametrization of two-doublet potentials 



In Section ^J] we wrote the two Higgs doublet potential as Eq. [|. Here we write 



two common forms of the most general two Higgs doublet potential. Firstly we 
write 

V(0i,0 2 ) = "7,?</>l0i + m 2 2 (t)\(j)2 + m\ 2 (j)\(l)2 + mi2020i 

+tri\<h<l>\<h + rA<h<&<h., (78) 

where the only complex parameters are the m\ 2 , £5, £&, and £7. This potential has 
14 independent parameters. Imposing the discrete symmetry 0i — > 0i, 02 ~ > ~<p2 
on dimension four terms will force £§ — if — 0, and we have a potential with ten 
independent parameters. 

Writing the same potential as 

Vfafa) = (A 1 + A 3 )(0 t 1 1 -^) 2 + (A 2 + A 3 )(0 2 2 -|) 2 

2 2 
+2A 3 (0 t 1 1 -^)(0 2 > 2 -^) 
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(79) 

where all the parameters are real, we again have a potential with 14 independent 
parameters. Imposing <pi — ► (f>i, (j) 2 — > —<f) 2 on dimension four terms we force 
four of these parameters H\ = \i 2 = /i 3 = /i 4 = 0, and we have a ten parameter 
potential. 

The advantage of writing the potential as Eq. |7U| is that the three of the 
parameters of the potential are £, v±, and -u 2 , and that the zero of the potential 



43 



is 



Vg 

V2 



10) 



where <p\ = 0, and <p 2 = £■ 

The relations between the parameters of Eq. 



and those of Eq. ^ are 
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-(A 2 + A 3 )t| - X 3 v 2 - ^-viv 2 cos£ - -£-v x v 2 sin£, 



ft. 



A, 



A 7 



^1^2 cos £ —V]V 2 sin £ — 



/ii 2 ft 2 



A 5 . . 
- — - WlU2 sin^ 

?i = Ai + Ag, 

h = A2 + A3, 
^3 = 2A3 + A4, 

A5 + A6 . 



4 
Ay 
4 



2' 



■Ui-u 2 cos£ 



ft. 



ft 



2- 



^(As - A 6 
~(ft -ift), 
7: (ft - «ft)- 



iA 7 ), 



81) 
82) 

83) 

84) 

85) 
86) 
87) 

88) 

89) 

90) 

91) 
92) 



We are free to redefine the fields Q of Eqs. [T8| and 79|. Rewriting Eq. [T9| with 



b a e liPa gives 



7(0i,02) = (A 1 + A 3 )(0l0i-^) 2 + (A 2 + A 3 )(0 2 02 -^) 



2 \2 



l>? n , ,+ , V 2 



+2A 3 (0 t 1 1 -^)(0 2 2 -^ 



+ X 4 



1^10202 



Re 2 (0l0 2 ) 



lm 2 (0 t 1 2 ) 



+ (A + + X i)(Re(0l02) - ^) 2 + (A+ - X i)lm(0 t 1 2 ) i 



+ X2 (Re(0 t 1 2 ) - ^)lm(0 t 1 2 



+ft(0|01 



V 



+ft(0^02-y)(Re(0 t 1 2 ) 



V1V2, 
2 ' 

V X V 2 , 



+ ft(0 t 1 1 -^)lm(0 t 1 2 



(93) 



9 / +ft(0^0 2 - ^)Im(ri0 2 ), 
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and we now have a potential which is a function of 13 parameters, one less than 
both Eqs. [78| and [79|. Where these new parmeters are in terms of those of Eq. |79| 



A+ 


2 

= ^(A 5 + A 6 ), 


(94) 


A_ 


— ^ 5 ~ 6 '' 


(95) 


Xi 


A 7 

= — sin 2^ + A_ cos 2^, 


(96) 


X2 


A 7 

= — cos 2^ — A_ sin 2^, 


(97) 


fix 


= /ii cos£ + /i 2 sin^, 


(98) 


ih 


= — nx sin£ + /i 2 cos£, 


(99) 


fiz 


= /i 3 cos£ + /i 4 sin£, 


(100) 




= — /i 3 sin£ + /i 4 cos£. 


(101) 






(102) 



On imposing the discrete symmetry <pi — + <fti, <p 2 — ► —<p2 on dimension four 
terms fii — fii — fii — fii — 0, and we have a potential which is a function of 
nine parameters, again one less than the potentials of Eqs. [75| and [79] with the 
same symmetry imposed. This nine parameter potential is Eq. ^| of section [Tl 
and is the potential we use throughout. 



B Extrema of the potential 

Extrema of the potential given in Eq. |9] occur at solutions to the four independent 
equations 

5V ^ .0, (103) 



(104) 



where i = 1,2, 3, 4, and Xi are the xi, x 2 , y%, and z 2 of 

, _ v\_ r o 1 = V2_ r z 2 

1 ~~ y/2 [ xi I ' 2 a/2 L x 2 + m 

A general bounded function of four variables with quartic and quadratic terms 
only can have up to 2 4 minima. 

The trivially found solutions to Eq. |103| are x\ = ±1, x 2 = ±1, y 2 = z 2 = 
(i.e. Eq. [LT]), and x x = x 2 = y 2 = z 2 
to find analytically was 



Xl 

2 i 2 i 2 i 

x 2 + y 2 + z 2 - 1 



x 2 



0. The only other solution we were able 

= o, 



A, 



(A 2 + A 3 )tan 2 /T 



(A+ + Xi) 



2/2, 



(105) 
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these describe a circle with one zero eigenvalue, and potential energy 



2„,2 



V 



ViV. 



l u 2 



AiA 2 + (Ai + A 2 )A S 
(A 2 + A 3 )tan 2 /3 



+ A+ + xi 



(106) 



9L 40, and 42, and is 



which may be less than zero for a potential obeying Eqs. 
a zero of the other terms of the static energy functional ^6|. 

To find numerically the global minimum, we implemented two methods. Firstly, 
using the Maple extremisation routine extrema, we looked for an extremum of 
V(Xi) with negative energy somewhere in the chosen region of parameter space. 
As the vacuum in our parametrisation has zero energy, this meant it was not the 
global minimum. We used this solution as an initial configuration for a simple re- 
laxation algorithm, which is equivalent to setting £ of the Newton method (Eq. 
[5^| ) to unity. We then scanned though parameter space relaxing to the global 
minimum at every point. 

Our second method was to use an initial configuration of X{ = 0, find the 
eigenvalues of the configuration, and add a perturbation in the direction of the 
most eigenfunction with the most negative eigenvalue. We then used the relax- 
ation routine on this configuration. We did this for each point in parameter space, 
reinitialising to X{ = at each point. 



C Static energy functional 

On substituting the ansatz of Eqs. |43|-f45| into the Lagrangian |2| we obtain the 
static energy functional of Eq. [K]. Here we give the form of K , K\, Vq, V\, and 
V2 for the C conserving ansatz and for the C and P violating ansatz. 

In the absence of C violation F a = a a and G a = c a , and we have the usual 
ansatz of Ratra and Yaffee (42| where K\ — V\ — and Kq, Vq, and V2 are 



Vn 



16M,1 



K + K 



G 



K D - — 

" 2r 2 

1 



1 



'22 



'22 



a a r + c a r + a 



(107) 
(108) 



K G 
^0 



2r 2 L4r 2 

1 -(al + cl)(a 2 + p 2 + 2) 



V2{3, 



- cl) - y/2. 



2 



(109) 



(Ax + A 3 ) (a 2 + c 2 -4cos 2 /?)' 

+ (A 2 + A 3 ) (a 2 + c 2 -4sin 2 /3) 2 

+2A 3 (a 2 + cl - 4 cos 2 0) (a 2 + c\ - 4 sin 2 (3) 

+A 4 (a x c 2 - a 2 ci) 2 
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+ (A + + xi) ( a i a 2 + cic 2 - 4 cos (3 sin f3) 2 



(110) 



Vo 



-A 4 + A+ - xi) (oic 2 - a 2 ci)" 



(111) 



This ansatz will maintain spherical symmetry if V 2 = 0. The condition V 2 = is 
met if A4 = A+— xi, or equivalently if M H ± = Ma- In cases where M H ± 7^ Ma, the 
spherical symmetry of a field configuration will still be maintained if aic 2 = a 2 ci, 
as the V 2 terms vanish from the energy density The ordinary sphaleron comes 
into this class of configurations since c\ = c 2 = 0. However, it is still important 
to include this term as it affects the form of £ used in Eq. 56 to calculate the 
curvature eigenvalues. 

In the presence of C violation b a and d a are no longer zero and Kq, K\, Vq, 
Vi, and V 2 are 



Vo 



Vi 



2r 2 



16M, 2 , 



K D 
-"■0 


+ K G 


•"■0 


1 - 
~ 2^ . 


K G 
^0 


1 " 
- ^ _ 



^(« 2 + f 2 -2) 2 



'2^2 



'2 



+ i («i + ^ + 4 + a(« 2 + /5 2 + 2) 



j'2 
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(113) 



v/2/3 



( a o + b 2 a - c 2 - d 2 a ) - V2a(a a c a + b a d c 



(a' a d' a - b' a c' a )r 2 + - (a a d a + b a c a )(a 2 + [3 2 - 2) 

(Ai + A 3 ) (a\ + bl + c\ + d\-A cos 2 p) 2 
+ (A 2 + A 3 ) (a 2 2 + b 2 2 + c 2 + d\- 4 sin 2 0) 2 

+2A 3 {a\ + b\ + c 2 + d\ - 4 cos 2 /?) (a 2 , + fr 2 , + c\ + c/ 2 , - 4 sin 2 0) 
+A 4 ((aic 2 - a 2 ci + 61 d 2 - 6 2 c?i) 2 

+ (aid 2 + a 2 d\ - b\_c 2 - b 2 cif - 4 {a\di - b\Ci) (a 2 d 2 - b 2 c 2 )) 
+ (A + + xi) (a-io-2 + hb 2 + c\c 2 + d\d 2 - 4 cos (3 sin (3) 2 
+(^+ - Xi) («i&2 - «2&1 + cid 2 - c 2 di) 2 

+2x 2 {aia 2 + b\b 2 + cic 2 + did 2 - 4 cos/? sin /3) (ai& 2 - a 2 &i + cid 2 - c 2 di) 

4(Ai + A 3 ) (a 2 + 6 2 + c 2 + d\ - 4 cos 2 0) (Mi - 61 Ci) 

+4(A 2 + A 3 ) (a 2 + 6 2 + c 2 , + d 2 , - 4 sin 2 /3) (a 2 rf 2 - 6 2 c 2 ) 
+4A 3 ( (a 2 + 6? + c\ + rf 2 - 4 cos 2 /3) (a 2 d 2 - & 2 c 2 ) 

+ (a 2 + 6 2 + c 2 + d 2 - 4 sin 2 0) (a x di - 6 x ci)) 
+2(A + + Xi) (o-i«2 + M2 + cic 2 + rfid 2 - 4 cos /? sin 0) (a x d 2 + a 2 di - b x c 2 



(114) 
(115) 
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— 2(A+ - Xi) (M2 - Mi + cid 2 - c 2 di) (aic 2 - a 2 ci + b ± d 2 - Mi) 

+ 2 X2 [(oia 2 + &i&2 + cic 2 + M2 - 4 008/98111/9) (aic 2 - a 2 Ci + M2 - b 2 d 1 ) 

— (M2 — Mi + M2 — Mi) (M2 + Mi — &1C2 — b 2 ci)} , 

4(Ai + A 3 ) (Mi - 6i Cl ) 2 

+4(A 2 + A 3 ) (M 2 - &2C2) 2 

+ 8A3 (Mi - &lCl) (M2 - &2C2) 

-A 4 ((aic 2 - a 2 ci + M2 - Mi) 2 
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-2% 2 (M2 + a 2 d x - b x c 2 - b 2 Ci) (a x c 2 - a 2 c\ + b\d 2 - b 2 di) J . 



D Numerical scheme 



(117) 



(118) 



To implement the scheme numerically, we discretise the n fields into N values 
Jai in the range < r < R. The values at the boundaries fAQ and /a(at-i) are 
determined by the boundary conditions in a way which we specify below. Hence 
E" is a n(N — 2) x n(7V — 2) matrix, and 5/ and £ ' are n(iV — 2) column vectors. 

To increase the accuracy of the solution while minimising the number of points 
N we use a rescaled co-ordinate s, where 



In C 



In 



1 + fir 
1 +r 



C 



1 + 

1 + i? ' 



(119) 



Here, M ma:r is the maximum of [Mh, Mh, Ma, Mh±], and for M max = M w we used 
M max = 1.01 x M w . We took R = 20M" 1 and used TV = 51 points throughout. 
It is also convenient to define two new functions X(s), Y(s) through 

1 1 (C s - /i) 2 



X(s) 



ds 
dr 



ln|C|(/i-l) C s 
1 (C s -/i)(C s + /i) 



H-l C s 

The first derivative of the energy £' may be split into Higgs and gauge parts 



(120) 
(121) 



S'c 



.(Yr 2 + 2r) 



dfn 
ds 



2 d 2 f H 1 d 
Xr , - + 



ds 2 Xdf H 



+ V + -V 2 



_ Y dfc_ x lA 
ds ds 2 



a + 1 d 



Xdf G 



Kg + V + ^V 2 



I 

3 



, (122) 
(123) 
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We use symmetric second-order accurate differencing for the derivatives, and so 



f 



(y. r 2 _|_ 2r )— Hi+1 ~ ^ m ~ 1 > v „2 (fm+i — 2fm + fm-i) 



1 d 



Hi 



2h s 



Kg + V 0i + '-v. 



XiTi 



hi 



2i 



(124) 



c Gi 



v (/Gi+1 - fGi-l) v (/Gi+1 - 2/Gi-l + /Gi-l) 
■Ji : ^-i 



1 d 



Xidf t 



Gi 



h a 



Kg + V 0i + '-V 2i 



hi 



3 



(125) 



where the index % — 1, (iV — 2), runs over the rescaled co-ordinate s, excluding 
the first and last points, and h s = (N — l) -1 , is the separation between each 
adjacent rescaled co-ordinate. We did not use (fm+2 ~ 2 fm + fm-2) I '(2h s ) 2 for 
the second order derivative, as this would have produced two systems independent 
in derivative terms, one seeing the even points and one seeing the odd points. 
The matrix £ is a block tridiagonal n(N — 2) x n(N — 2) matrix of the form 



^i-l,i-2 



D T,i-i 



i+l,i+2 



where each of these boxes are n x n matrices, and there are (N — 2) x (N — 2) 
such boxes. The only non zero terms are the D'^, D^, and D^ i+1 . Note that 
and D^ i+1 are themselves diagonal, with entries 



D~ - 1 Y - 1 X 



'2/1, * 



(126) 
(127) 



for the two gauge fields, and 



D 



^(2. + ^)-^, 



(128) 
(129) 
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for the remaining Higgs fields. If we write 

u Ai,Bi — - 

then D®f er are diagonal in A, B with 



r^Oder 
u ALBi 



u Ai,Bii 



j-\0der _ 

i,i 


2 Y 


(gauge fields) 


-ptOder 

i,i 




(higgs fields). 



The non-diagonal elements are symmetric in A, B with 



D 



Omat 
Ai,Bi 



d 2 



dfAidf, 



Bi 



Kg + Voi + =rV z 



2i 



(130) 

(131) 
(132) 

(133) 



We have to be careful about the form of £ ' at the top left corner of the ma- 
trix, corresponding to the i — 1 point, affecting the D xx , and the D X2 terms. 
Also the bottom right corner, corresponding to the i = (N — 2) point, affect- 
ing the -D^_ 2iV _ 3 , and the D%_ 2N _ 2 since these must implement the boundary 
conditions. 




D 



N-3,N-4 



D 



N-2,N-3 



D 



N-3,N-2 



u N-2,N-2 



Because for the sphaleron the boundary conditions at the origin are never 
updated, D X1 , and Df 2 f° r the sphaleron are as Eq. |130| . For the RWS, and 
bisphalerons at the origin we use for the gauge fields 







fc\i=0 — /g|z=1) 



/(3lr=0 

from this we are able to calculate 

*|i=o = arctan(-a//?)| i=0 . 



(134) 



(135) 
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For the Higgs fields we use Tables |3] and | to give 



c a |j = o — a Q |j=o tan 0Q,|j =o , (136) 
d a \i=o = b a \ i=0 tan Q a \i=o, (137) 

where the Q a \i=o are calculated from \l/|j = o of Eq. |135| , and using Tables § and 
§ according to whether we are looking for the bisphalerons or RW sphalerons. 
Further imposing smoothness of (fr^p at the origin gives boundary conditions 





i=0 


= a a 


i=x cos 2 6 Q 


\i=0 + c a 


i=1 sin0 Q 


i=0 cos a 


j=0; 


(138) 


b a 


i=0 


= b a 


i= i cos 2 Q 


i=0 + 


i= i sin a 


i =0 COS Q 


i=0j 


(139) 


Qv 


i=0 


Cq, 


i= i sin 2 Q 


i=0 + 0>a 


i= lCOS0 Q 


i=o sin Q 


i=0) 


(140) 


d a 


i=0 


da 


i= i sin 2 a 


i=0 + b a 


i= l COS Q 


i=o sin 9 Q 


8=0- 


(141) 



To update the origin after each Newton Raphson iteration we use Eqs. |134 



138| - |141[ , We also use these to give us the form of D® 1 and D± 2 when looking for 



the bisphalerons or RW sphalerons. We did this by first writing, for a a - 

(\r 2 , r, \ dda I v 2 d Cta i 

— {Yr +2r)——\i—Xr 



ds IX ds 2 
1 

2hs 



-(Yir 2 + 2n)— — (a a \ 2 - a a \i cos 2 Q |i - c a \i cos6 a |i sinGali) 



-X x r\— z(a a \2 - 2a a \i + a a U cos 2 a |i + c a U cos 6 Q |i sin a |i), (142) 

with the equivalent expression for the other Higgs fields; and using Eq. |134| , for 
the gauge fields, we write 

- Y t i' - x ^ = - y ^(^i= - /°w - - feu- < 143 > 

We then, after functional differentiation of Eqs. |142| and |143| , get a form of l 
and Df 2 that sees the boundary conditions. 

The a throughout are zero if we are looking for sphaleron solutions, and are 
determined from either Tables ||| or § with Eq. |135| according to whether we are 
looking for bisphalerons or RWS. 

We now turn to the boundary conditions at infinity. The last point is never 
updated since this boundary does not evolve, and D° N __ 2N _ 2 is as Eqs. [130| - |133 . 



We did not use fo\r-^oo — (/g|jv— i — fG\N—2)/hs — as the boundary condition 
since rescaling the radial co-ordinate to allow greater accuracy at the origin re- 
duces the number of points at large distances. This meant that the form of the 
first and second derivative were not very accurate at the last few points. 

The form of £ of Eqs. [T24] and was not affected by the boundary con- 
ditions. Because £ is only defined for i = 1, (N — 2) and first and second 
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derivatives at i = 1, and i = N — 2 are obtained from the already updated fields 
/a|o and /a|jv-i- 

Also recalling that £ of Eqs. 56 and 57 used in the evaluation of the curvature 
eigenvalues is functionally differentiated with repect to fc and rfn, and not fa 
and fa- The form of D\ 1 and Df 2 for evaluating the curvature eigenvalues is for 
the Higgs fields components as Eqs. |130| and |133| since 5(rfn)\o = 0. We again 
use Eq. |134| for the gauge fields. 

To find solutions other than the original sphaleron we first find the sphaleron 
and determine the curvature eigenvalues and eigenfunctions of the configuration. 
If there is more than one negative curvature eigenvalue, we succesively add a 
fraction of the eigenfunction of the second (or third) negative eigenvalue to the 
sphaleron field configuration, measuring the energy at each step. If we chose this 
fraction small enough (typically between 0.01 and 0.1) the energy at each step 
will decreases until it reaches a minimum. When the energy after a step is larger 
than the energy measured after the previous step, we multiply the fraction by 
—0.1 and continue until the fraction is — 10~ 9 times its original value. 

This configuration is then used as the initial configuration for the Newton 
Raphson minimisation routine to find the RW sphalerons (or bisphalerons) . 

Sliding down the most negative eigenfunction of a sphaleron configuration 
reaches the vacuum. Sliding down the second most negative eigenfunction reaches 
the lowest energy branch of sphaleron like solutions, a third negative eigenfunction 
will reach the second lowest energy branch and so on. In this way we were able 
to find bisphalerons and RW sphalerons of the theory. 

We use BLAS fortran subroutines dgbco and dgbsl to solve for 5f a of Eq. |56] 
and subroutine dgeev to evaluate the curvature eigenvalues and eigenfunctions. 
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